Decoherence of intermolecular entanglement in exchange-coupled nanomagnets 
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We theoretically investigate the hyper fine- induced decoherence in a pair of spin-cluster qubits, 
consisting of two exchange-coupled heterometallic wheels. We identify two distinct regimes in the de- 
coherence of intermolecular entanglement and show that this can be substantially recovered through 
dynamical decoupling. Different chemical elements and physical processes are responsible for the 
decoherence of the singlet-triplet superposition, resulting in a wider tunability of its decoherence 
time. 
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Molecular nanomagnets (MNs) represent a rich class of 
spin clusters, whose properties can be widely tailored by 
chemical synthesis [2 J. Weak interactions between well 
defined molecular units have also been demonstrated [3] 
and can be made tunable by introducing different inter- 
molecular linkers [4] . This has recently allowed to demon- 
strate equilibrium-state entanglement between pairs of 
heterometallic wheels pQ. Quantum mechanical features, 
such as free [5j [6] and forced [7H9] oscillations, have also 
emerged in the coherent dynamics of MNs. The tunable 
coupling between nanomagnets, combined with the co- 
herent manipulation of their spin state by pulsed EPR, 




FIG. 1: (color online) (a) Chemical structure of the (CrrNi)2 
dimer 1 , as determined by X-ray diffraction: Cr (purple), Ni 
(green), F (yellow), O (orange), C (black), N (cyan); the H 
nuclei have been omitted for the sake of clarity, (b) Schemat- 
ics of the spin system and environment. The state of the two 
electron-spin rings (faint blue), is mapped onto that of two 
exchange- coupled 1/2 spins, Sa and Sb (blue). The electron 
spins interact with a bath of N n = 312 nuclear spins (red). 



will possibly enable the experimental investigation and 
the controlled generation of intermolecular coherence and 
entanglement. Besides being of fundamental interest, 
these issues are pivotal to the implementation of quan- 
tum information processing with MNs [T0lfT3] . 



The main sources of electron spin decoherence in MNs 
are represented by intermolecular dipole-dipole inter- 
actions [14], spin-phonon [15] and hyperfine coupling 
[16] [17]. While the former two contributions can be 
limited by dilution of the molecular ensemble and low 
temperatures [5], the suppression of hyperfine-induced 
decoherence might require more sophisticated dynami- 
cal decoupling techniques [18]. In spite of its funda- 
mental and applicative interest, still little is known on 
intermolecular coherence in coupled nanomagnets and 
on its robustness with respect to hyperfine interactions. 
Here we address this problem in the case of two spin- 
cluster qubits, each consisting of a Cr^Ni heterometal- 
lic wheel (Fig. fl[a)) [19 . Being only the ground-state 
doublet of each wheel involved in the low-temperature 
dynamics, the electron spin state of the (Cr7Ni)2 dimer 
can be mapped onto that of two exchange-coupled spins 
Sa — Sb — 1/2 (Fig. [IJb)). Essentially, decoherence is 
due to the dependence of the slow nuclear bath dynam- 
ics on the state of the electron spins, and to the resulting 
correlations between electron and nuclear spins. The fea- 
tures that characterize the decoherence of different linear 
superpositions can however remarkably differ from one 
another [20]. We indeed find similar differences between 
the decoherence of two prototypical linear superpositions 
(|3>i) + | $2)) /V% with \$i) the electron-spin eigenstates 
of the dimer. The first linear superposition, with compo- 
nents |$i) = I ^aAb) and |$i) = | ^aAb), corresponds 
to an entangled state of the two Cr^Ni molecules. Here, 
electron-nuclear correlations are shown to arise mainly 
from the opposite chemical shifts of the H nuclei induced 
by |3>i) and |3>2)- I n the second case, the components 
are instead |$ 1>2 ) = (| ^aAb) ± I H>aAb))/V2- Both 
the singlet and triplet are characterized by a vanishing 
expectation value of the electron spin projection, and 
therefore induce no chemical shift of the nuclei. We find 
that, as a consequence, the electron- nuclear correlations 
result from processes that are second-order in the hyper- 
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FIG. 2: (color online) Simulated time evolution of the en- 
tanglement between Sa and Ss, as quantified by the con- 
currence, Cab(Pc ) (black circles, left axis). The entangle- 
ment between the electron and nuclear spins, for a randomly 
generated state of the nuclear bath |I), is quantified by the 
von Neumann entropy (red circles, right axis). Figure inset: 
calculated time-dependence of n (black circles) and fitting 
curves f(t) = exp{ — (t/rd) a } (gray lines), with a = 3.31 and 
Td = 0.70 /xs and a = 1.46 and ra = 0.75 /xs. 



fine interaction (i.e. electron-spin mediated transitions 
between pairs of nuclei) and selectively involve the F nu- 
clei, localized close to the electron spins. 

The model — In the (CrrNi) 2 dimer, the spin system 
and environment consist of N e = 16 electron spins and 
N n = 312 nuclear spins, respectively (Fig. [lib)). Each 
CryNi wheel is modeled by a ring of eight electron spins - 
one for each of the magnetic ions - with S1-7 = SQr = 3/2 
and s$ = SNi = 1. The two lowest eigenstates of each 
single-ring Hamiltonian form a Kramers doublet: | ff^ = 
\S X = 1/2, M x = 1/2) and | ^ = \S X = 1/2, M x = -1/2), 
where \ — A B. At zero magnetic field this doublet is 
separated from the first excited quadruplet by an energy 
gap of A c± 20 K pQ. At temperatures T«A, the state 
and the couplings of each ring can thus be mapped onto 
those of an effective 1/2 spin. This also applies to the 
intermolecular exchange induced by the bipyridine linker 
in the (Cr7Ni) 2 dimer (Fig. [ljb)) PQ, and to the hyperfine 
couplings to the nuclear spins [21] . In the presence of an 
applied magnetic field the effective two-spin Hamiltonian 
reads: 



1-Lab = Jab$>a • S# + g[iB 2_^ ^ 



B, 



(1) 



=A,B 



with Jab <^ A. The spin bath consists of N n = 312 
nuclear spins (6 F and 306 H, with J F = J H = 1/2), 
whose positions are derived from X-ray crystallography. 
The nuclear spins I& interact with each other and with 
the electron spins via dipole-dipole interactions, which 
are accounted for within the point-dipole approximation. 
Energy relaxation of the electron spin state, aris- 
ing from transitions between different eigenstates |S, M) 
(with S = Sa + S#), is strongly inhibited, being the dif- 



ferences between the eigenvalues of Hab are much larger 
than both the excitation energies of the nuclear spins 
and of the hyperfine couplings. In such pure-dephasing 
regime, the electron-spin decoherence is induced by the 
dependence of the nuclear-spin dynamics on the electron 
spin state, and by the resulting electron-nuclear correla- 
tions. In particular, an initial state 

|*o> = |* c o)®|*no) = ^(|5i,Mi) + |52,M 2 ))(»|I) (2) 

(where |\P e o) and |^no) = |I) = (S^Jifc) are the electron 
and nuclear components, respectively) evolves into 

1 



l*t>= 
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|5i,Mi)®|Ii(t))+e^W|S r 2,M 2 )®|l2(t)) • (3) 



Here, |Ii(t)) and |l2 (^)) are the states of the nuclear bath 
conditioned upon the electron spins being in the | Si , Mi ) 
or |S2, M 2 ) states, respectively. The reduced density ma- 
trix of the electron spins is thus given by: 



p e (t) = Tr n {\9 t )(%\} = ± J2 \S k ,M k )(S k ,M k \ 

fc=l,2 

\S 2 ,M 2 )(S u M 1 \e i ^ +h.c], (4) 



n(t) 



where r\ = | <(Ii 1 12 ) | and <pi = arg{(Ii|l2)} result from 
the dynamics of the nuclear bath. 

This is derived from the effective Hamiltonian [2TJ [22] : 



Heff =(T e z ®Hi e) +l e ®Hil\ 



(5) 



being \ty t ) = exp{-zH e ff£/ft}|^o). Here, the 

electron and nuclear-spin operators, are given by: 
a% = |5i,Mi)(5i,Mi| - |S 2 ,M 2 )(S 2 ,M 2 |, X e = 
|5i,Mi)(5i,Mi| + |5 2 ,M 2 )(5 2 ,M 2 |,and 
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^i a) = E4 a) 4* + E \BShii? + c%)i+ir , (6) 
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with a = e, i. The intrinsic nuclear Hamiltonian Wa' in- 
cludes the nuclear Zeeman and the secular (i.e. Zeeman- 
energy conserving) terms of the dipole-dipole interactions 
between nuclei. The extrinsic Hamiltonian Hn accounts 
for the Overhauser effect (couplings A^) and for the in- 
teractions between nuclei mediated by virtual excitations 
of the electron-spins (#£ and C^ ). 

Bell-state decoherence — In order to investigate the 
decoherence of intermolecular entanglement, we consider 
the case where the components of |^ e o) (Eq. [5| are 
Si = Mi = 1 and S 2 = — M 2 = 1, which corre- 
sponds to the spins S^ and S# being in the Bell state 
|*eo) = 7^ (lfU,lh?) + \UAb)) • The simulated build 
up of correlation between the states of electron and nu- 
clear spins and the resulting loss of entanglement between 
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FIG. 3: (color online) Time evolution of n, averaged over 
N — 100 randomly generated initial states |I), in the ab- 
sence or in the presence of a 7r-pulse applied to the electron 
spins at t = r (solid-gray and dotted curves, respectively; 
different colors correspond to different values of r). The val- 
ues (n(t = 2t)) for different r (red squares) are fitted by 
/(£) = exp{ — (t/rd) a }j with a = 1.8 and Td — 2.9 /xs. 



S^ and S# is reported in Fig. [2] The entanglement be- 
tween two qubits is quantified by the concurrence [23] , 
which, for the density matrix p e (Eq. [3]), reduces to: 



CAB(pe) =rj. 



(7) 



The von Neumann entropy [23] can be used to quantify 
entanglement between electron and nuclear spins. Its 
expression, for \Sfr t ) (Eq. |3|, turns out to be: 
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As p e evolves from the initial linear superposition to a 
statistical mixture of |^a?1Tb) and | J|a 5 JJ-s), Cab (black 
circles) and S en (red circles) tend to and 1, respec- 
tively. This clearly shows the mutually exclusive nature 
of electron- nuclear and intermolecular entanglement. Be- 
sides, two distinct regimes can be clearly identified in the 
entanglement decoherence: for times t < 1/is, rj(t) fol- 
lows an exponential decay f(t) = exp{ — (t/r^)"}, with 
a = 3.31 and r^ = 0.70 /is; for t > 1 /is, instead, a = 1.46 
and Td = 0.75 /is (figure inset). 

For temperatures much larger than the nuclear Zee- 
man splitting, the density matrix of the nuclear bath 
resembles the maximally disordered mixture of states 
|I>: Pno = Ei|I)(I|/2^. While n(t) is largely in- 
dependent on the initial state |I) of the nuclear bath, 
this is not the case for the phase </>i, which depends on 
the Overhauser field induced by the nuclei (see Eq. [6]): 

B^ = - Y,k=i A-k Ik/(PB9z) (g z is the effective g factor 
of each CryNi ring) . The static component of the Over- 
hauser field and the resulting inhomogeneous broadening 
produce a fast dephasing, that can however be eliminated 
by spin-echo techniques [5]. The slow dynamics of the 



nuclear bath however also results in a fluctuating com- 
ponent of the field, with a potentially irreversible con- 
tribution to the phase <j>i(t) and to the electron-spin de- 
coherence. The simulated dynamics however shows that 
such contribution affects only marginally entanglement 
decoherence in the (Cr7Ni)2 dimer: 



CAB(pf") = (rie^)^( ri ) 



i4>i 



(9) 



where < \{rie i(t)1 )\ - (n) < 0.012 for t < 10 /is, and 
the average is performed on a set of randomly generated 
initial states |I). 

In order to gain further insight into the decoherence 
process and underlying dynamics, we isolate the role 
played by the specific chemical elements and the different 
physical mechanisms. To this aim, we consider the differ- 
ence between Kne^ 1 )! resulting from the full dynamics 
and the same quantity obtained by neglecting part of the 
terms in H e ff- As detailed in Table [I] the contribution 
of the H (F) nuclei to the decoherence of intermolecular 
entanglement is predominant (negligible). Amongst the 
hyperfine coupling mechanisms, the effect of the chemi- 
cal shift (i.e. the renormalizations of the nuclear Zeeman 
energies induced by the electron spins) largely prevails 
over the electron-spin mediated interactions between the 
nuclei. The relative importance of these contributions 
changes drastically in the case of the singlet-triplet deco- 
herence (see below). 

The correlations between electron and nuclear spins de- 
scribed above can be partially reversed by dynamical de- 
coupling schemes [18 J. In Fig. [3] we show for example the 
effect on entanglement decoherence of a flip of the elec- 
tron spins S x= a,b at a time t = r, that can be induced by 
an EPR 7r-pulse. This ideally transforms the two compo- 
nents of |* e0 ) into one another (\1[aAb) < — > \$aAb)), 
and thus c e z — > —&% in H e ff at t — r. Such inversion 
of the extrinsic part of the Hamiltonian tends to undo 
the electron- nuclear correlations; (n) correspondingly in- 
crease for t > r, with a maximum localized around 
t ~ \[2t (colored lines) [24]. The above mentioned tt- 
pulse is used within the Hahn-echo sequence, in order 
to cancel the effect of inhomogeneous broadening in free- 
induction decay experiments and to induce an observable 
spin echo at time t = 2r. At that time, electron- nuclear 
correlations have built up again (red squares), determin- 
ing values of (rj) which are lower than those obtained in 
the absence of the tt— pulse (brown curve). More com- 
plex dynamical decoupling schemes are thus required in 
order to achieve simultaneously the removal of inhomoge- 
neous broadening and electron-nuclear entanglement and 
to make the latter effect experimentally observable [25 . 

Singlet-triplet decoherence — In order to complement 
the results on the entangled state decoherence, we con- 
sider the case of a fundamentally different linear combi- 
nation, namely that between the singlet and triplet states 
(Si = 0, S 2 = 1 and Mi = M 2 = in Eq. [2j. The re- 
duced density matrix, derived from Eq. [3] now reads: 
p e (t) = |[(l+r)| UAB)(UAB\ + (l-r)\ UAb)Aa 
Ab I ] (with r = (rie^ 1 )), and is always factorizable into 
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FIG. 4: (color online) Singlet-triplet decoherence as a function 
of the intermolecular exchange Jab the Zeeman energy Ez — 
gHBB z (Eq. 111). The function r = (ne^ 1 ) has been averaged 
on N — 128 randomly generated initial states |I). 
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TABLE I: Main contributions to the decoherence of the Bell 
state (BS) and of the singlet-triplet superposition (ST). The 
key quantity is Sr = max t |(rj) — (n)|, where rf is computed 
by keeping in H e f£ (Eq. [5} only the chemical elements and 
the terms specified in the central columns; the considered 
time interval is < t < 10 jjls. 



the states of S^ and S#: Cab(Pc) = 0. 

The features of the singlet-triplet decoherence are fun- 
damentally different from those emerged in the case of 
the Bell state. As shown in Fig. |4j the time-evolution of 
r is strongly dependent on the intermolecular exchange 
Jab (Eq. pi). In particular, the decoherence rate is on 
average much smaller than that of the Bell state, and de- 
creases for increasing energy difference between the sin- 
glet state and each of the triplet states. An analogous 



dependence (not shown here) is found with respect to the 
magnetic field. These features are fully consistent with 
the identification of transitions between F nuclei, medi- 
ated by virtual transitions of the electron-spin state, as 
the dominant contribution to the overall decoherence (see 
Table]!]). The dynamics of the H nuclei, though present, is 
largely independent on the electron-spin state, and there- 
fore doesn't result in the built up of the electron nuclear 
correlations. As anticipated in the introduction, such cor- 
relations arise in the case of the entangled state from the 
fact that the two components (| ilA^ilB) and | ^a^b}) 
give rise to opposite chemical shifts. Instead, the chemi- 
cal shift induced by the singlet and M = triplet states 
vanishes, being (S,0\sf z \S,0) = 0. The dependence of 
the nuclear spin dynamics on the electron state therefore 
relies essentially on second-order processes and on the 
contribution of the F nuclei, that are localized closer to 
the electron spins of the dimer (Fig. []Ja)). 

Conclusions — While a number of additional aspects 
need to be considered in view of their possible relevance 
in experiments (for example, those related to the hyper- 
fine interaction with the nuclei in the solvent, or to the 
intialization and manipulation of the electron spin state 
of the dimer), some fundamental features have emerged 
from the presented simulations. In particular, we have 
identified the clear occurrence of two regimes in the build 
up of electron- nuclear correlations, and of the resulting 
decoherence of intermolecular entangled. The compar- 
ison with the case of singlet-triplet decoherence shows 
that the decoherence in the low-energy subspace of the 
dimer can be induced by distinct chemical elements and 
physical processes, depending on the components of the 
initial linear combination. Finally, as demonstrated in 
other spin systems, dynamical decoupling might allow to 
disentangle electron and nuclear spins, thus leading to a 
substantial recovery of the intermolecular coherence. 
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